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Abstract 

Any nonsingular function of spin j matrices always reduces to a matrix polynomial of order 2 j. The 
challenge is to find a convenient form for the coefficients of the matrix polynomial. The theory of biorthogonal 
systems is a useful framework to meet this challenge. Central factorial numbers play a key role in the 
theoretical development. Explicit polynomial coefficients for rotations expressed either as exponentials or as 
rational Cayley transforms are considered here. Structural features of the results are discussed and compared, 
and large j limits of the coefficients are examined. 


1 Introduction 

As a consequence of the Cayley-Hamilton theorem mm any nonsingular function [3] of an n x n matrix always 
reduces to a polynomial of order n — 1 in powers of the matrix, since the nth and higher powers of the matrix can 
be reduced to a linear combination of lower powers (see [5] and earlier literature cited therein). The challenge 
for anyone wishing to take advantage of this fact is to find a convenient form for the coefficients of the matrix 
polynomial. 

In particular, a nonsingular matrix function of any one component, ft • J , for a unitary, irreducible angular 
momentum representation of spin j, always reduces to a matrix polynomial of order 2 j. Explicit polynomial 
coefficients for rotations expressed either as exponentials or as rational Cayley transforms mm are considered 
here, to augment some recent studies euzhbi. Structural features of the results are discussed and compared, 
including the behavior of the coefficients for large j. 

The theory of biorthogonal systems musim is a useful framework to find the coefficients of the matrix 
polynomial, even in the most general case. For analytic spin functions of ft • J, central factorial numbers nam 
play a key role in the construction of the biorthogonal system. 

In prior work, Curtright, Fairlie, and Zachos (CFZ) obtained explicit and intuitive results [B] expressing the 
standard, exponential rotation matrix for any quantized angular momentum j as a polynomial of order 2 j in the 
corresponding (2 j + 1) x (2 j + 1) spin matrices ft • J that generate rotations about axis ft. While many previous 
studies of this or closely related problems can be found in the literature nai mama on Dans!, none of these other 
studies succeeded to find such simple, compact expressions for the coefficients in the spin matrix polynomial, as 
elementary functions of the rotation angle, as those obtained by CFZ. For each angle-dependent coefficient in 
the polynomial, the explicit formula found by CFZ involves nothing more complicated than a truncated series 
expansion for a power of the arcsin function. Although a detailed proof of the CFZ result is not exhibited in 
j6), the essential ingredients needed to provide such a proof are in that paper, and indeed, the details of two 
elementary derivations were subsequently given in [7 . 

More recently, Van Kortryk [8] discovered the corresponding polynomial result for the Cayley rational form 
of an irreducible, unitary SU (2) representation. At first glance the Cayley transform for a spin representation 
would seem to follow immediately from the CFZ result just by changing variables from 9 to a function 9 (a), where 
a parameterizes the transform. For any given numerical value of ft • J this would be so, obviously, but it is not 
obviously so for matrix-valued ft ■ J. Nevertheless, as it turns out, the result for the Cayley transform is actually 
much simpler than the CFZ result for the exponential. For the Cayley transform the explicit coefficients involve 
nothing more complicated than truncated series expansions of already finite polynomials in a real parameter a. 
Remarkably, this is true for the Cayley transform of any matrix, not just spin matrices union nil- 

in this paper I compare the CFZ and Van Kortryk results, and I provide a more complete discussion of the 
latter, including a detailed derivation of the polynomial coefficients for the SU( 2) Cayley transforms. I show 
explicitly how the CFZ results are connected to Cayley transforms by Laplace transformations. 
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2 Methodology 

The theory that underlies construction of spin matrix polynomials is discussed in several places in the literature 
(see m [3 and references therein). Here I condense and summarize the key results. 

2.1 A fundamental identity 

The fundamental identity for spin j G {0, 1, |, • ■■ } is [7] 

2 j 

(h • jf J+1 = - ^2 1 ( 2 + 2j, 1 + m) x (n • J) m , (1) 

m—0 

where the coefficients t ( n , k) are central factorial numbers [T^JES]- (For details, see Appendices A and B.) This 
identity ensures that all powers of (n • J ) higher than 2 j can be reduced to linear combinations of lower powers. 
This is a specific illustration of the Cayley-Hamilton theorem, applied to the monomial on the LHS of ©• 


2.2 Vandermondeum 

The Vandermonde matrix [2^ for spin j is 


V\j} = 


1 2 j (2 jf 

1 2j — 2 (2j — 2) 2 


(2 jf ~ 

( 2 j - 2 ) 2j 


- 1 -2 j (-2 jf 


(-2 jf J 


( 2 ) 


while its inverse [23] is given by the “Fin du monde” matrix elements 




2 1 "* N k (l,j) 

(2J + 1-0! (i-1)! 


( 3 ) 


where k,l G {1, 2, • • • , 2j + 1}, and where the numerator factors Nk ( l,j ) are given by nested sums subject to 
some exclusions. 

! 1 for k = 2 j + 1 

X) (j + 1 - mi) ■ • • (j + 1 - m 2 j+i-k) for 1 < k < 2j + 1 ' 

l<m 1 <m 2 <---<m 2 j+i-k< 2 j+l 
mi,m 2 ,--- ,m 2j+ i-k^l 

( 4 ) 

These are unfamiliar polynomials [25| in l and j. For k = 2j + 1 the numerator is always 1, as it says above, 
and for k = 2j it reduces to a single sum giving a result only first order in both j and l. Thus 


N 2j +i (l,j ) = 1 , N 2j {l, j) = 1 — l + j . 


( 5 ) 


But for other values of k these numerator polynomials quickly get out of hand. For example, for k = 2j — 1 the 
resulting double sum yields a polynomial 3rd order in j and 2nd order in l. 

N 2j -1 (1,3) = (1 - if + 2 (1 - l)j + l (1 -j) (2 j - 1 ) j . ( 6 ) 

For reasonable values of j — up to 100, say - it is easiest to just use numerical machine computation to 
obtain V~ l [j] directly from V [j], with no need of the analytic closed-form expressions for the Nk (l,j)- 


2.3 Biorthogonal matrices 

It is useful to have in hand the dual matrices which are trace orthonormalized with respect to powers of the spin 
matrix, S = 2 h ■ J. Without loss of generality, choose S = 2 J 3 , since any other choice for h merely requires 
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selecting a different basis to diagonalize the spin matrix, thereby obtaining the same eigenvalues as 2 J 3 . 
the powers are 


(2 j) m 0 • • • 0 0 

0 ( 2 j — 2 ) m ■■■ 0 0 



0 0 ••• (—2j + 2) m 0 

0 0 ••• 0 (-2 j) m 


Thus 


( 7 ) 


It is now straightforward to construct orthonormalized dual matrices T n such that 


8 n ,m = Trace (T n S m ) , n,m = 0,1, • • • , 2j . ( 8 ) 

The T n may also be chosen to be diagonal (2 j + 1) x (2 j + 1) matrices in the basis that diagonalizes S. In fact, 

for any spin j the required entries on the diagonal of T n are just the entries in the (n + l)st row of the inverted 

Vandermonde matrix, V~ l [j]. (Note that here, unlike the conventions in [ 6 ], both rows and columns of the 
Vandermonde matrix and its inverse are indexed as 1, 2, • • • , 2j + 1.) That is, 

(Tn-i) kk = ( V - 1 [j]) n k , n, k = 1 , • • • , 2j + 1 . (9) 

This result follows immediately from the fact that the diagonal entries for S m are just the entries in the corre¬ 
sponding column (i.e. the (m + l)st column) of the Vandermonde matrix, V [j]. 

Explicit examples of spin matrix powers and their duals, and the corresponding Vandermonde matrix inverse 
V -1 [ 7 ], are given in Appendix C for j = 1/2, 1, 3/2, and 2. 


2.4 Projecting coefficients 

Suppose / is a nonsingular matrix function for spin j. In a basis where S = 2 J3, the matrix / ( S ) is diagonal 
and so is its reduction to a polynomial in S: 


f (2j) 0 

0 / ( 2 / — 2 ) 


0 0 
0 0 


0 0 
0 0 


/ (- 2 J + 2 ) 0 

0 / (— 2 j) _ 

2 j 

= E /™ = E 

m =0 m —0 


f(S) 


( 2 j r 

0 


0 

0 


0 

(2 J - 2) m 


0 

0 


Trace orthogonality of the dual matrices then gives the coefficients as 


0 0 
0 0 


(—2 j + 2 ) m 0 

0 (-2 j) m 


fm = Trace (T m f (S)) 

Assembling the coefficients into a column and using © gives 


/o 

fl 


E.l + i 1 (P- 1 [J]) lfc /(2J + 2-2fc) 

Efcii 1 (P _1 [J]) 2j fc / (2J + 2 — 2 k) 

= v - 1 [j] 

fm 

f (2j — 2) 

hj- 1 
hi \ 


Efcti 1 (V-'\j]) 2 . k f( 2 j + 2 - 2 k) 

_ Efcfi 1 (V~ 1 \fl) aj+lik f( 2 j + 2 - 2 k) 


/ (2 — 2/) 
/(—2j) 


( 10 ) 


( 11 ) 


So the spin matrix polynomial coefficients for / (S') are determined by straightforward matrix multiplication 
using the inverse Vandermonde matrix. 
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In terms of the explicit form in ([3]), for k = 1, 2, • • • ,2j + 1, 


-l 


2j+l 

/*-i = E (v 

i=i 

2 j+l 


>kl 


f(2j + 2-2l) = 2 1 ~ k 


i- fc /(2J + 2-2J) 


Z=1 


( 2 j + l-l)\ (l- 1 )! 


= 2 


l-fe 


E 

1=1 


(-i) 


k—l—2j 


f(2j + 2-2l) 


(2j + l-l)\ (l-l)\ 


E 


(j + 1 - mi) ■ ■ • (j + 1 - m 2 j+i-k) ■ ( 12 ) 


l<m 1 <---<m 2 j+i-k<2j+l 


The coefficients of the highest powers appearing in the spin matrix polynomial are often worked out readily from 
the general closed form expressions for the Nk ( l,j ), as given in (0 and 0. For the highest three powers in the 
polynomial, 


f 9 - 2 j ST' (- 1) 1 1 f ( 2 j + 2 - 21) 

123 (2J + 1-0I (i-1)! ’ 

f __ 9 i- 2 j ST' (- 1) 1 1 (1 — l+j) f(2j + 2 — 2l) 

/2J_1 E (2J + 1 - 0 ! a - 1 )! 

o2 _ 2 . (- 1 ) 1 "' ((1 -l) 2 + 2(1-1) j+ 1(1-j ) (2j - 1 ) j) / (2 j + 2-21) 

h >- 2 ~ 2 J E (2j + 1 - 0! a - 1)! 


(13) 

(14) 


(15) 


etc. But again, for reasonable values of j it is probably easiest to just use numerical machine computation to 
first obtain V[j] directly from V [j]. and then perform the matrix multiplication in ( 1111 ) . 

Indeed, by taking several low values of j and their explicit V~ x [j]. Van Kortryk [ 8 ] was able to deduce the 
general formula for the coefficients in the Cayley transform polynomial, which he then confirmed for specific 
higher values of j. I provide a proof of the formula in Section 4 of this paper. 


2.5 Lagrange-Sylvester expansions 

I note in passing another useful method, albeit equivalent to that involving the inverse Vandermonde matrix. 

Reduction of nonsingular matrix functions for spin j to polynomials can be efficiently carried out in specific 
cases using Lagrange-Sylvester expansions [2B]. Consider functions of an N x N diagonalizable matrix M with 
non-degenerate eigenvalues Ai, i = l,--- ,N. On the span of the eigenvectors, there is an obviously correct 
Lagrange formula, as extended to diagonalizable matrices by Sylvester, 

N 

/(M) = ^/(Ai)Pi, (16) 

i =1 

where the projection operators — the so-called Frobenius covariants — are given by products, 


N 

n>, = n 


3 =1 


M - A j 
A i A j 


(17) 


From expanding each such product it is evident that any /(M) reduces to a polynomial of order N — 1 in powers 
of M, 

N -1 

/ (M) = £ C m [/] , (18) 

m =0 

and the function-dependent coefficients can be expressed in terms of the eigenvalues of M by expanding the 
projection operators ED as polynomials in M. 

CFZ used the methods described above to obtain the coefficients for the spin matrix polynomial reduction 
of the exponential form for rotations for several low values of j. They then re-expressed the results in a form 
that surprisingly suggested an immediate generalization, after which they developed the theory to confirm their 
ansatz, thereby raising it to the level of a theorem with explicit proofs [7]. I review and illustrate the CFZ 
formula in the next section. 
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3 Rotations as exponentials 


For comparison to the results of the next section, we first recall the CFZ formula for a rotation through an 
angle 0 about an axis ft, valid for any spin j £ {0, i, 1, |, • • • }. The formula reduces the manifestly nonsingular 
exponential function of any spin matrix to an explicit polynomial: 


2 j 


exp(i 6 h ■ J) = 4 ? 1 (0) (2 i fi ■ J) k , 


k —0 


where the angle-dependent coefficients of the various spin matrix powers are given by 


(9) = sin fc (9/2) (cos (9/2)) e ^' k ^ Trunc 

\j~k/2\ 


f arcsinVx 


(vi - V 


x=sin 2 (S/2) 


(19) 


( 20 ) 


Here, |_- • • J is the integer-valued floor function and Trunc [/ (x)] is the nth-order Taylor polynomial truncation 

n 

for any f (x) admitting a power series representation: 


oo n 

f®) =52 fmX m , Trunc [/ (x)] = 52 / m x m . (21) 

z—/ n z—/ 

m =0 m—0 

In addition, e ( j , k) is a binary-valued function of 2 j — k that distinguishes even and odd integers: e (j, k) = 0 
for even 2 j — k, and e ( j , k) = 1 for odd 2 j — k. 

The simplest two nontrivial cases of (fTTFl) are well-known: j = 1/2 and j = 1. Explicitly, 

exp(i 6 n ■ J)\ j=1/2 = A 1 q /2] (9) J 2x2 + 2sh4 l i 1/21 (9) (n ■ J ) 2x2 

= cos (0/2) I 2 x 2 + 2* sin (0/2) (n • J ) 2x2 , (22) 


exp (i 0 n • J )| i=1 = 4 1] (0) 1 3 x 3 + 2 ZH ' 11 (0) (n • J ) 3x3 - 2 kl' 1] (0) (n • J)Ls 

= -f 3 x 3 +*sin( 0 )(n-J ) 3x3 +(cos 0-1) (n-J ) 3x3 . (23) 

The former involves the Pauli matrices, J = er/2, while the latter is sometimes known as the Euler-Rodrigues 
formula. Several other explicit cases may be found in PH and [ 6 ]. 

In practice, for finite j of reasonable size, the truncations needed to evaluate (1201) are easily obtained as a 
matter of course by machine computation, for example by using either Maple® or Mathematica®. Nevertheless, 
it is interesting and useful for a systematic analysis that Taylor series for powers of cyclometric functions can be 
expressed in terms of t (to, n), the so-called central factorial numbers of the first kind [12;, 13]. Thus for \z\ < 1 
and non-negative integer n (cf. Theorem (4.1.2) in |13|). 


(arcsin (z)) n 


m=n 


(24) 


Note that the coefficients in these Taylor series are all non-negative. In general, the values of t (to, n ) are defined 
by and obtained from simple polynomials, as described in Appendix A. 

Incorporating (l24l) into the expression for the coefficients (l20l) gives 


1 .1 o m 

A® (9) = ^y —r \t (to, k) I sin " 1 (0/2) for even 2 j — k. (25) 

2 fc ml 

m=k 

As firmly established in nam, the remaining coefficients in (1191) may then be obtained from 

4-1 W = | ^ 4 1 ( 0 ) for odd 2j-k + l. (26) 

Also, as observed in [ 6 ], the results (1201) display the limit j —> oo for fixed A: in a beautifully intuitive way. In 
that limit, the truncation is lifted to obtain trigonometrical series for the periodicized 9 k monomials. But even 
as j —> oo, integer j (bosonic) and semi-integer j (fermionic) coefficients are clearly distinguished by a relative 
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sign flip for 9 £ [tt, 37t] mod (47r). This is evident upon plotting the first few coefficients for very large spins. For 
example, Aq\ ... 5 ( 9 ) are plotted here for j = 69 (darker curves, in blue) and j = 137/2 (lighter curves, in red). 








As stated in (1261) . the slope of each red (fermionic) curve plotted on the right is given exactly by the red curve 
plotted to its immediate left. Similarly, the slope of any blue (bosonic) curve on the left is given by the blue 
curve to its right, but in the row above. 
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4 Rotations as Cayley transforms 


The Cayley rational form of a unitary SU (2) group element is also nonsingular for any irreducible representation 
with spin j. Therefore it can also be reduced to a spin matrix polynomial of order 2 j. Thus 


1 + 2 ia n • J 
1 — 2 ia ft • J 


^fc 1 ( 2 * h • J ) fe 


k—0 


(27) 


Here a is a real parameter, to avoid singularities for any j € {0, 1, §, 2, • • • }. The coefficients 2lj^ (a) 

are to be determined as functions of a. These coefficients can be obtained by rewriting the geometric series 
1/ (1 — 2 ia ft • J) for spin j as a polynomial in n • J. Thus define the auxiliary polynomial 


1 


2 j 


77 = E ®* 1 (<*) ( 2 * «• J ) fc ■ 


1 — 2 * 0 ! n • J 


k—O 


Then clearly 


01 b] — — 1 01 bl — 


The simplest two nontrivial cases are again j = 1/2 and j = 1. Explicitly, 


1 


1 — 2io n • J 


= S ' 1721 (o) / 2x2 + 2 * (o) (n • J) 


3 = 1/2 


[1/2] 


1 


l + O 2 


(/ 2x2 + 2 io (n • J) 2x2 ) , 


(28) 


(29) 


(30) 


1 + 2*o n • J 


1 — 2*o h • J 


= *4 1/2] (o) J 2x2 + 2 * 2l/ /z| (o) (n • J) 


i=i/2 


[ 1 / 2 ] 


2x2 


1 + 0 ' 


((1 - o 2 )/ 2 x 2 +4io (n-J) 2x2 ) , 


(31) 


= ^i , 11 (a) J 3x3 + 2 * S™ (o) (n • J ) 3x3 - 4 ®' 1] (o) (n • J) 2 x3 
= Y 7-77 ((! + 4 « 2 ) J 3X3 + 2*0 (n ■ J) 3x3 - 4o 2 (n • J) 3x3 ) > ( 32 ) 


1 + 2*o n • J 


1 — 2*o h • J 


= a ' 11 ( a ) I 3x3 + 2 *; a‘ iJ (o) (h • j ) 3x3 - 4 a‘ iJ (o) (n • j) 


.[ 1 ] 


[ 1 ] 


*3x3 


J=1 


1 


1 + 4o 2 


((1 + 4a 2 ) J 3x3 + 4*0 (n • J ) 3x3 - 8 a 2 (n • J) 3x3 ) 


(33) 


Yet again, the first two of these polynomials involve the Pauli matrices, J = <r/2, while the last is related to the 
Cayley transform of the Euler-Rodrigues formula. Appendix D repeats these two examples and also gives the 
Cayley transforms for j = 3/2, 2, 5/2, and 3. 

The coefficients in the geometric series expansion for arbitrary integer or semi-integer j follow directly from 
the methods in 13 0- The result is [ 8 ] 


(«) = -j—71—— -77 Trunc [det (1-2 ia h ■ J)\ , 

det (1 - 2*o n • J) 2 j-k 

where the truncation is in powers of a, and where the determinant for spin j is 

2 j L/+1/2J 

det (1 — 2*o n • J) = (1 — 2*o (j — m)) = ^1 + 4o 2 ( j + 1 — r*) 2 ^ . 

m—0 n= 1 

These results are readily checked for small values of j upon using explicit matrices, say n • J = J 3 . 
detailed proof of ( 1 M 1 ) in the following subsection. 


(34) 


(35) 
I provide a 
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But first, note for integer spins it is always true that ©{/^ (a) = 1, and that the higher coefficients are paired, 


© 


L j] 

2fc+2 


(a) = a © 2 fc+i ( a ) f° r 0 < fc < j — 1 and integer j , 


while all the coefficients are paired for semi-integer spins. In this case 

© 2 fc+i ( a ) = Q ® 2 fc ( a ) f° r 0 < k < j — i and semi-integer j . 


(36) 


(37) 


So, the coefficients in the spin matrix polynomials for bosonic (integer j) and fermionic (semi-integer j ) Cayley 
transforms are easily distinguished, both quantitatively and qualitatively. The coefficients are paired as indicated 
because, for either bosonic or fermionic spins, only even powers of a are produced by the determinant factors in 


(1311) . Hence the parities: 21^ (—a) = (—l) K 2l^ J (a) and ©J? 1 (—a) = (—1) K ©^ (a). 

Moreover, for any allowed j all the coefficients of a 2k in (1351) are positive. As a consequence, for any allowed 
j the 2l[^ (a) and ©[^ (a) coefficients have no singularities for real a. In particular, for any j the highest two 
coefficients reduce to 

i ®0 («) = — ' ' 1 


k oi bl 


ill] 


\ k oi bl 


73T ®S-x (a) =- 


(38) 


a 2 i 2j v ' a 2 -! -1 2j 1 v ’ det (1 — 2 ia h • J) 

And as j increases, the domain where 1/det(l — 2 ia n • J) achieves significant values collapses towards the 
origin, at which point it always has unit value. Here is a plot for j = 1/2, 1, 3/2, 2, 5/2, & 3. 



1/ det (1 — 2 ia n • J) plotted for j = 4,1, |, 2, and 3, for upper to lower curves, respectively. 

Two additional comments are warranted [ 8 ]. First, the determinants in (1341) are essentially generating 
functions of the central factorial numbers t(m,n) (see Appendix A), a fact previously exploited in [ 6 ) 7] and 
already emphasized for the Cayley transforms in [ 8 ]. For either integer or semi-integer j, 

B—fe/2j 

Trunc [det (1 — 2 ia n • J)j = 4 m a 2m |f ( 2 j + 2, 2 j + 2 — 2m)| , (39) 

2j—k ^^ 

m =0 

with the full determinant obtained for k = 0 if j is integer. For semi-integer j . the same expression is true 
for truncations, only now the full determinant includes an additional term 4 J+ 2 a 2j+1 1 1 (2 j + 2,1)|. Thus the 
general formula for the determinant, for either integer or semi-integer j, is 

b'+jJ 

det (1-2 iaft-J)= ^ 4 k a 2k \t { 2 j + 2, 2 j + 2 - 2 k)\ . (40) 

fc=0 

Second, as j —> oo for any fixed k the truncation in (13-H) is lifted — but with some subtleties to be discussed 
below — to obtain for small a, lim^oo ©^ (a) ~ a k . But in contrast to the periodicized 0-monomials found 
in [ 6 ], the large j behavior here does not make the periodicity of rotations manifest. To exhibit periodicity even 
for finite values of j, 6 must be expressed, on a case-by-case basis, as cyclometric functions of a, and then the 
branch structure of those cyclometric functions must be invoked ( 8 ). 








4.1 A derivation for general matrix resolvents 

In this subsection I obtain the coefficients in the Cayley transform spin matrix polynomial by deriving the 
expansion coefficients for the resolvent of a general matrix. This is a well-known result in resolvent theory (e.g. 
see [ 2 Q]) and has been established many times in the literature (e.g. see mmy While there is hardly any 
need to plow again a field that has been tilled as much as this one, I include here a variant of the standard proof 
just to make the discussion self-contained. (Also see Appendix D.) 

For any finite N x TV matrix M consider the resolvent written as a matrix polynomial, 


1 


JV-l 


1 — aM 


= '^2 r m (a) M r 


m —0 


This is much simpler than the exponential case in that one need not differentiate, e.g. to find 

d ( 1 \ M 

(1 - aM) ■ 


da V 1 — olM 


1 — aM 


(41) 


(42) 


so as to obtain an expression linear in M multiplying the original polynomial, and subsequently a set of first-order 
equations for the coefficients r m . Rather, it suffices here just to multiply through by 1 — aM. Thus 


JV-l 


1 = (1 - aM) = (1 - aM) £ M m r m . 


m =0 

On the other hand, the Cayley-Hamilton theorem gives 

N N -1 

v M m d m = o, m n = -— v M m d m , 

dN ^0 

where the polynomial coefficients d n are defined by the determinant 

N 

det (1 — aM) = ^ Oi m d m , 

m =0 


with do = 1, djv = (—1) W detM, etc. (see [27]). Therefore C51) becomes 


(43) 


(44) 


(45) 


1 = ro + 


ar N - i 
dN 


N -1 

d 0 +J2M m 

m =1 



arm -1 


ar jy-i 

dN 



(46) 


If the powers of M m for 0 < m < TV — 1 are all linearly independent, as is the case for the spin matrices h ■ J, 
then 

ro + Qr j V ~ 1 do = 1 and r m — ar m -1 + 01 f -1 d m = 0 for 1 < m < N — 1 . (47) 

UN dN 

It is now straightforward to show the solution to these recursion relations is given by 

a n 

r n = -——-— Trunc [det (1 — aM)) . (48) 

det (1 — aM) N-i-n 1 v n v ' 

For M = 2 in • J with TV = 2j + 1, the r n are essentially just renamed * 8 ^. Thus the result (1341) is proven. 


4.2 Additional exact results and asymptotic behavior for large j 

The determinants can also be written as Pochhammer symbols: x^ = x (x + 1) ■ ■ ■ (x + n — 1) = r ^”- > . For 
example, for integer j, 

n (1 - 2ia O' “ TO )) = n 0 - 2ian ) = (-2 ia) J+1 r , (49) 

m =0 n=1 V 2 iot) 

’■I i r (i 4-1 4- j- j 

(1 — 2 ia (j — m)) = ]^[ (1 + 2 ian) = (2ia) J+1 — , 1 . 2za , (50) 

m=j -(-1 n =1 \2iaJ 
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which leads to a nice form for the determinant expressed in terms of analytic functions of both j and a. Again, 
for integer j , 

\i+i r 0’ + 1 ~ 2 k) F U + 1 + 2 k) 


det (1 — 2ia h • J) = (4a 2 ) 


r(-sfe) 


r(^) 


= — (2 a) 2 ^ +1 sinh 
7 r 




( . , 1 \ 

( 2 ^ 

X 

v + + 2ia) 


upon using T (+—) T (— 515 ') = 2na/ sinh On the other hand, for semi-integer j, 


det (1 — 2 ia n • J) = — (2a) 


2j+l 


cosh 


(a 


r '’ + 1 + 2 S 


(51) 


(52) 


In the limit j —> 00 , these are the well-known infinite product representations of sinh and cosh. That is to say, 


lim 
j ^°° (4a 2 ) 


det (1 — 2 ia fi ■ J) 


b'+lj 


(r(i + ^)) 2 


- sinh (t^;) for integer j 

■ cosh (^) for semi-integer j 


(53) 


This shows that there is a distinguishable difference between Cayley transforms for bosonic and fermionic spins, 
even as j —> 00 , much as there were distinguishable differences in the behavior of the coefficients for the expo¬ 
nential m, also as j —> 00 , as illustrated in the Figures of Section 3. 

Now consider the truncations for k = 1 or 2, say for integer spin. This choice of coefficients removes from 
det (1 — 2 ia n ■ J) the highest power of a, namely the a 2 - 7 term. But this term is easily computed: 


n= 1 

Therefore for integer j 


J J 

^4a 2 (j + 1 — n) 2 ^ = (4a 2 fc 2 ) = 4 J a 27 (j !) 2 = 4- 7 a 2 - 7 |T (j + 1 )| 2 . 


(54) 


fc= 1 


Trunc [det (1 — 2 ia n ■ J )] = Trunc [det (1 — 2 ia n • J)] = det (1 — 2 ia n ■ J) — A^a 2 ^ (jiy 

2j-l 2j-2 


= 4 j a 2j ( — sinh x T (j + 1 + 

I 7 r V2a/ \ 2 1 a 


- I r (i + 1)1 


(55) 


The corresponding coefficients in the Cayley rational form are 

a 


(a) = 


det (1 — 2 ia h ■ J) 
= a [ 1 - 


—— ^det (1 — 2 ia n • J) — 4Pa 2 ^ 
4 J a 2 ^' |T ( j + 1)| 2 


(j'0 2 ) 


i ( 2 a) J sinh (£) x |r(j+l + ^)| 


(56) 


Similarly for (a) = a (a). The exact results for these two coefficients, for all integer spins, are then 
given by 

_ 7 T l r (j + !)| 2 

^ 'nr ■ 1 t v 1 — 1 — - 


1 (a) = T *2 ] (a) = 1 - 


a 2 asinh(^) x |r(i + 1 + 5 ^) 

For fixed a, the large j limits of these follow immediately from 


(57) 


lim 

j^-oo 


i r (j + i)r 


= 1 . 


(58) 

l r (j + 1+ 2^)1 

In fact, the limit is just 1 to all orders in 4-. This should be obvious, but in case it is not, as a check the reader 
may consider the following expansion: 

\ 2 

1 - Psi (1 ,j + 1 )z 2 + — {6 (Psi (1, j + l )) 2 - Psi (3 ,j + 1)) 


( r U + 1)) 


r (j + 1 + z) r (j + 1 - z) 


( 59 ) 


^ (Psi (5, j + 1) - 30 Psi (3, j + 1) Psi (1, j + 1) + 60 (Psi (1, j + l)) 3 ) 


360 

1 


70 (Psi (3 ,j + 1))" + 840 (Psi (1, j + l )) 4 - Psi (7, j + 1) 

20 160 V +56 Psi (5, j + 1) Psi (1, j + 1) - 840 Psi (3, j + 1) (Psi (1, j + l )) 2 

+ 0(z w ) , 
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where Psi (n, z) is the nth derivative of the digamma function Psi(z) (also known as PolyGamma in Mathemat- 
ica ®). Except for the O ( z°) term, all of the O ( z k ) coefficients in this last expansion vanish as j —> oo. 
Therefore it follows for integer j that 


T («) = 3 <*f 0 


1 


1 - 


j->oo 2asinh(^) 


(60) 


So, as a approaches zero, the RHS approaches unity faster than any power of a due to the essential singularity 
in the exponential, 

(61) 


7r 7T 

— e 2 RT 


2asinh(^) <*-»o |a 
This is the precise meaning of the statement in B) that 

lim * 8^1 ( a ) ~ a 1,2 . 

’ o-t 0 


j->oo 


(62) 


In fact, the asymptotic form is rapidly approached as j is increased. 
A similar calculation, again for integer j, gives 


A ' («) = "4 1 (“)=!- 


^ l r (j + i)r 


a 


1 - 


1 + 


j-*oo 2asinh(^) V 24a 2 ) 
For the case at hand, 


2 asinh(^) x |r (j + 1 + ^) 
_.2 


Xi? ( 1+ 24^ ( 7r2 - 6Psi ( 1 » 1 +i))) ( 63 ) 

(64) 


Psi(l,l+i) = C(2)-Efc2 • 


fc=l 


(65) 


Here are some other facts that were used to get exact expressions as well as asymptotic behaviors for the 


* 8 ^ 1 4 (a) coefficients: 


f[(^ 2 ) - 


1=0 


r (j +1 - \Xi) r (j +1 + -/X) 

r(-vX)r(vX) 


1 1 ( 2 j + 2 , 2 )| — — (z + l 2 ) 


1=0 


= (r (j + 1 )) 2 


z=0 


1 d 


2 J 


|t( 2 j + 2,4)| - g ^2 II ( z + l 2 ) 


1=0 


z=0 


= g (T (j +1 )Y n 2 - (r (j + i)y Psi ( 1,1 + j) , 


( 66 ) 

(67) 

( 68 ) 


Trunc 2 j-i [det (1 — 2ix h ■ J )] Trunc 2 j -2 [det (1 — 2 ix n • J)] ^ A^x 2 ^ (T (j + l )) 2 


det (1 — 2 ix h ■ J) 


= 1 - 


det (1 — 2 ix h ■ J) 

( r (j + l )) 2 

231 H^)l r (j + i + i)| 


det (1 — 2 ix n • J ) 


(69) 


Trunc 2 j -3 [det (1 — 2ix n ■ J )] Trunc 2 j -4 [det (1 — 2ix h • J )] ^ A^x 2 ^ (T (j + 1 )) 


det (1 — 2 ix h • J) det (1 — 2 ix h • J) det (1 — 2 ix h • J ) 

y-igpj-i ^1 (r (j + l )) 2 7 T 2 - (T (j + l )) 2 Psi (1,1 + j)) 


det (1 — 2 ix n ■ J) 


= 1 - 


k l r (i + 1 )! 2 


2 asinh(^) x |r (j + 1 + 


r+ ( 1 + i( l2 ‘ 6Psi(1 ’ 1+j) )) 


(70) 
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I plot a few examples to illustrate these features, beginning with 18 


[j] 

1 


(a) and *8^ (a). 



03^ (a) /a = *8|f J (a) /a 2 plotted for j = 1, 2, & 8 (lower to upper black curves) and the j —> oo limit 

(uppermost curve, in red). 


ibl 


A closer look near a = 0 is provided by the next graph. 



03^ (a) /a = © 2 J (a) /a 2 for j = 1, 2, 




& 8 (lower to upper black curves) and the j —> oo limit (uppermost 
curve, in red). 
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I also plot the corresponding results for *83^ (a) and *84^ (a). 



183 ^ (a) /a 3 = © 4 1 (a) /a 4 for j = 2 , 


di] 


4, & 12 (lower to upper black curves) and the j 
curve, in red). 


00 limit (uppermost 


Again, a closer look near a — 0 is given in the next graph. 



183 ^ (a) jo? = 2 > 4 J (a) /a 4 for j = 2 , 




4, & 12 (lower to upper black curves) and the j 
curve, in red). 


00 limit (uppermost 


Similar behavior is exhibited by the higher coefficients. 
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These plots should be compared to the large j behavior of the coefficients for the exponential, given in (BUI) . 
As the Figures in Section 3 attest, the asymptotic j —> oo behavior is not closely approached by the ( 9 ) for 
fixed k until rather large values of j are considered, and for some 6 not even then. This is especially significant 
for 6 near ±7r, where Hindoo ( 0) has a discontinuity in 9 for every other value of k. By contrast, for fixed k 
as j —> oo the 21 (a) converge uniformly and monotonically to (a) on any finite interval in a, and indeed, 

the large j limit of the coefficients is closely approached for moderate values of j. These last points are clearly 
displayed in plots of the relative error between 21 (a) and (a) for a few values of j. Further numerical 

studies suggest that, for any fixed k, the relative error 


Ai? 1 (a) = 


sir (a) - a!?' 1 (a) 


21 fc 1 (a) 


(71) 


is always > 0 and goes monotonically to zero as j —> oo, for all a. 
Here are some examples. 



The relative error (a) = (a) plotted for j = 1, 2, & 8, as upper to lower curves, respectively. 



The relative error A^ (a) = A^ (a) plotted for j = 1, 2, & 12, as upper to lower curves, respectively. 
Again, similar behavior is exhibited by the higher coefficients. 


14 
















Large j conjectures and proofs for all the coefficients 

Here are the first four requisite asymptotic j —> oo results, for integer j, as obtained by explicit calculation: 


asynqq (a) = 1 
asymp 2 (a) = 1 
asymp 3 (a) = 1 
asymp 4 (a) = 1 


1 


2a: 

7 r 

sinh 

7 r 

2a 


1 


2a 

7 r 

sinh 

7 r 

2a 


1 


2a 

sinh 

7 r 

7T 

2a 


1 


2a 

7 r 

sinh 

7T 

2a 


/ 1 

/ 7 r \ 

1 + — 

— 

1 3! 

V2a/ 

( 1 

/ 7T \ 


— 

1 3! 

V2a/ 

( 1 

/ 7T \ 


— 

1 ^3! 

V2a/ 


)*)• 


5! \2aJ 


1/7r \ 6 
7! V2 a) 


(72) 

(73) 

(74) 

(75) 


(Semi-integer j are different. See below.) Note that the numerators are just truncations of the series for the 
denominator, ^ sinh ^ ( 2 ra +i)! (la) 2 "- Hence the obvious conjecture for integer j is 

i /k—1 o \ / 7 r\l+ 2 fc 


asymp fc (a) = 1 —- 


sinh 


E 


/ 7r 


, 2 n 


(£) 


2 ol \n=0 


(2n+ 1)! V2a, 


(1 + 2k)\ sinh 


hypergeom I [ 1 ] 


2 a 


1 + k, — + k 


Given this result, the coefficients of the spin matrix powers in the reduction 

/fc-i 

ilbl _ ' ' 

fb;(2 n+l) 


16 a 2 J 
(76) 

are, as j —> oo for integer j, 


lim 

j—t oo a 


®“-‘ ( “> = ““ 7S W' 1 - 2a si ) h , E (STS (£) 


2 fc-l 


2n 


j ^oo a. 


for k > 1 . (77) 


2a \n=0 


As an immediate consequence, lirn^oo asymp fc (a) = 0 for all a / 0, with the limit approached monotonically 
from above, while at a = 0 , asymp fe ( 0 ) = 1 for all k. 

Also, to put it briefly, explicit calculation of the lowest few coefficients reveals the corresponding conjecture 
for semi-integer j: The asymptotic expressions may be obtained from the bosonic spin case by substituting 
cosh y-. Thus for the fermionic case, as j —> oo for semi-integer j , 


— sinh 

7 r la 


lim -jnr ( a ) = lim 2 Li ^ 2 fc+i 

J —>-00 J—>• OO ^AC-l-1 


<8 


[j] 


(a) = 1 — 


1 


E 


/ 7T 


cosh^ \ ( 2 n)! \2a 

2a \ n= o v > 


r), 


for k > 0 . 


(78) 


While plots of the various functions of a in m and (1781) are qualitatively similar, nevertheless, the bosonic and 
fermionic Cayley transforms are still easily distinguished, even as j —> oo. 

Proofs of these conjectures are straightforward using (13H) . (1391) . (l5.il) . and Proposition 2.7 in [15] . For 
example, for integer j, use the result (15i)l) to rewrite (13H) for the even-indexed coefficients as 


1 


k 


1 _L_ ( a \ = _ 

a 2k 2k det (1 — 2 ia h • J) 


E( 2a ) 


2(7+1 -0 


\t(2j + 2,2l)\ . 


The limit j —> oo now follows from the integer j result in (1531) , written as 


(2af(T(l+j)Y 


lim 


7 —ioo l det (1 — 2 ia n • J) J 2a sinh (^) 

and from the asymptotic form of |t (2 j + 2,2m)| as given in [13], Proposition 2.7 (xxvi), written as 


lim -* \t{2j + 2,2l)\ = 


1 


(21-1)! ( 2 a) 


2(/-i) 


(79) 


(80) 


(81) 


j^oo (r(i + j)) 

The net result is an expression whose numerator is a truncation of the infinite series for the denominator, as 
conjectured above. 

/k -1 


lim ( 1 - ® 2 fe («) ) = 

7—100 \ 


1 


sinh 


E 


1 


/ 7 r 


2 a \n —0 


( 2 n + l)! V 2 a 


r)- 


(82) 


The odd-indexed coefficients for the bosonic spin case are given by this same expression according to the ele¬ 
mentary pairing (1361) . Thus for integer j the conjectured j —> 00 form of the coefficients is established. 

A corresponding proof for semi-integer j follows from (1531) and Proposition 2.7 (xxvii) in [13] . 
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The result for semi-integer j is as conjectured in (1751) . Alternatively, it may be expressed as a hypergeometric 
function similar to that in (1761) . 


1 - 


1 


cosh 

la 



l 

(2 n)\ 



( _7T_\' 2 + 2 k 

V 2a) _ 

(2 + 2 fc)! cosh y- 


hypergeom 



3 

2 


/c, 2 T k 


16a 2 


(83) 


5 Relating exponentials and rational forms by variable changes 


To finish the comparison of rotations as exponentials and Cayley transforms, for arbitrary j, it is instructive to 
return to an issue briefly mentioned in the Introduction. Namely, why does a map from one form to the other 
not succeed through the use of a numerical-valued change of variable a (0) or its inverse 0(a)? After all, a 
comparison of the exponential and Cayley results for j = 1 / 2 , as well as for j = 1 , as given in (1551) . (1551) . (15T1) . 
and (1331) . would give the misleading impression that this can succeed. 

To see the difficulty, take a basis where h • J is diagonal, with eigenvalues M G {—j, — j + 1, • • • , j — 1, j}, 
and equate the exponential and Cayley forms acting on a specific eigenstate \M). The result of course is 


exp (i 9 M) 


1 + 2 ia M 
1 — 2*a M 


(84) 


This is easily solved either for a (6) or for the inverse map 9 (a): 


a ( 0 ) = tttt t an 
v ’ 2 M 



9{a) 


2 

M 


arctan( 2 Ma) . 


(85) 


So the relation is M-dependent, although in this particular instance it is independent of the sign of M. In 
general, nontrivial M-dependence would be a feature that would arise upon comparing any two functions of 
n • J when acting on | M). 

Now it just so happens that a unique map works for all the M -eigenstates for j = 1/2, and another unique 
map works for all the M-e igenstates states for j = 1. (Note that for M = 0 the relation (l84l) always reduces to 
1 = 1 no matter what the relation between a and 9.) But this is fortuitous and peculiar to those two spins due 
to \M\ having only one nontrivial value in those special cases. For all j > 3/2 the map must change with \M\. 
There is no consistent a 9 relation that works on all n • J eigenstates for j > 1. In a sense made precise by 
m , there is a “parameter shear” in the relation between a and 9 as the h ■ J spectrum is traversed. 

Similarly, if one were to equate the coefficients of any given power (n • J ) in the spin matrix polynomial 
reductions of the exponential and the Cayley transform, the result would be an a ^ 0 relation that depends on 
k, and this result would be inconsistent (except in the two special cases j = 1/2 and j = 1) with that obtained 
from equating the coefficients of some other power of n • J . 


6 Relating exponentials and rational forms by Laplace transforms 


Some standard results from resolvent theory provide a direct relation between (USD and (1551) . Consider the 
Laplace transform and its inverse, 

roo p'y+iT 

F (a) =£[/]= / e~ st f (t) dt , f(t)=C~ 1 [F} = — lim/ e st F(s)ds. ( 86 ) 

Jo 27r * T ^°° Jj-iT 

For any hermitian N x N matrix, a Laplace transform yields the resolvent from the exponential, and vice versa. 


1 


s — iM 


= £ [e ltM ] 


e itM = c -i 


s 


1 

— iM 


(87) 


This furnishes another route to derive the CFZ formula starting from the known matrix polynomial expression for 
the resolvent. Given the ease by which the polynomial coefficients for the resolvent are obtained, this provides 
perhaps the simplest proof of the CFZ result. 

For linearly independent powers M k , 0 < k < N, as is the case for spin matrices, the matrix polynomial 
expansion coefficients for the exponential and the resolvent are distinctly related order-by-order by the Laplace 
transform. 


JV-l 


exp ( itM ) = ^ A n (t) (iM) n , ^ B n (s) {iM) n , B n ( s ) = C [A n ] , A n (t) = C 1 [B n ] . ( 88 ) 


N—l 


71=0 


s — iM 


71=0 
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If the first k powers of M are not independent, for k < N — 1, correct relations between the exponential and the 
resolvent are still obtained by taking C or £ _1 , as in (1571) and (l88l) . although the final results may admit further 
simplification. Fortunately, it is not necessary to worry about that situation in the following. 

Rescaling variables produces a form of the Laplace transform that is more useful for the direct comparison 
of (P2U1) and (1M1) . Thus 


1 — 2 i a n • J 


e * exp (2* at n • J) dt 


> o 


1 c°° 

(a) = y J o e_t A k ] ( 2ai ) dt ■ 


(89) 

(90) 


It suffices here to consider only the direct transform. The inverse relation is automatic, given the large a 
behavior of the 58^ (a). From (1251) . along with (OBI) , the required calculation is 


1 

1$ (a) = — - e _t sin m (at) dt for even 2 j — k . 

J 0 


These are elementary integrals, of course, resulting in 


m/2 


I®{a)=a m H 


1 


( m + 1)/2 


Z=1 


1 + 4Z 2 a 2 


for even m , 1$ (a) = a m ]Q 


z=i 


1 + (21- 1)a 2 


(91) 


for odd m . (92) 


After combining these results with (1251) and (1351) . and some elementary manipulations of the central factorial 
numbers, the transform (1901) yields the expected result: 


®S? ] (a) = . E 4m « 2m I* & + 2, 2j + 2 - 2m) | . 

' m —n 


a 


U-k/2] 


det (1 — 2 ia h 


m =0 


For example, for spin j = 1/2, with the usual 2x2 matrices, 

exp (i 9 fi ■ J) = cos (0/2) I + 2i sin (0/2) (n • J) , 


/* oo /*oo 

/ e _t cos (at) dt = - -- , / e _t sin (at) dt = — 

J o 1 + oi 2 Jo 1 


+ a 1 


det (1 — 2 ia h ■ J) = 1 + a 2 , 


1 


1 


1 — 2 ia h ■ J 1 + a 2 
and for spin j = 1 , with 3x3 matrices, 


(I + 2ia fi • J) 


exp (i 0 h • J) = I + i sin (0) (h ■ J) + (cos 0 — 1) (n ■ J) 2 , 


/ e *dt = 1 
'0 jo 


j* OO poo 

, / e -t cos (2at) dt = -—^ , / e _t sin (2at) dt = 

Jo 1 + 4 « Jo 


2a 


1 + 4a 2 


det (1 — 2 ia n ■ J) = 1 + 4a 2 , 


1 


1 — 2 ia h ■ J 1 + 4a 2 


~^~2 ((l 4 a2 ) 4 + 2ia n • J— 4a 2 (n ■ J) 2 ^j . 


(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 
( 100 ) 
( 101 ) 
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7 Concluding remarks 

Polynomial reductions of spin matrix exponentials have long been recognized as important and useful in many 
different contexts m-m, em. Perhaps the results presented in Ena® and discussed further in this paper 
can help to facilitate these and other applications. 

Moreover, Cayley transforms have numerous applications, many of a practical nature. For example, consider 
[32] and [33]. In the first of these two papers, previous results on attitude representations were generalized using 
Cayley transforms, for application to guidance and control problems, while in the second paper, the theory of 
time stepping methods were developed in a systematic manner based on the Cayley transform, for application to 
discretised differential equations that describe evolution in Lie groups. Both of these papers stress that Cayley 
transform methods are easier to implement than the evaluation of matrix exponentials. 

This relative ease of implementation is borne out in the present paper by direct comparison of (1201) and (1341) , 
and more indirectly by the relative simplicity of the derivation for the Cayley transform coefficients, as given 
here, when compared to the proofs for the coefficients in the exponential case, as given in [7]. 

Acknowledgements I thank T S Van Kortryk and C K Zachos for discussions. I have also benefited from 
stimulating views of sea waves approaching the Malecon, La Habana, and I thank the organizers of STARS 2015 
for their efforts to make this possible. This work was supported in part by NSF Award PHY-1214521, and in 
part by a University of Miami Cooper Fellowship. 
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Appendix A: Central factorial numbers 


For historical reasons, central factorial numbers are defined as the coefficients in simple polynomials mm- 
They can be either positive or negative, but only their absolute values are needed for the coefficients of the spin 
matrix expansions in the main text. Moreover, t (even, even ) are integers, but t (odd, odd ) are not integers, and 
t (odd, even ) = 0 = t (even, odd). So the even and odd cases are best handled separately. 

By definition and as elementary consequences thereof (cf. http://oeis.org/Al82867), 


( x 2 — l 2 ) = ^ t ( 2 m, 2k) : 


,2k 


1=0 


k=l 


-k 1 d‘ 


k m— 1 


t ( 2 m, 2k) = (-1 ) m - fc |t ( 2 m, 2 fc)| = (-l) m “ fc ± ft (* + l 2 ) 


1=0 


z=0 


as well as (cf. http://oeis.org/A008956) 


n k- ( i + 


1=0 


^ t ( 2 m + 1, 2k -f 1 ); 


, 2 fc+l 


k =0 


t ( 2 m + 1 , 2k + 1 ) = (-l) m - fc 1 1 (2m + 1 , 2k + 1 )| = (-l) m “ fe A ^ J] ^ + (i + j 


z=0 


(Al) 

(A2) 


(A3) 

(A4) 


Appendix B: Proof of the fundamental identity 

Here I provide a demonstration that 

2 j 

(2 h ■ J) 2j+1 =-^2 2 1+2j ~ k xt( 2 + 2 j, 1 + fc) x (2n • J) k , (Bl) 

k —0 


where t (m, n) are the central factorial numbers, defined in Appendix A. In a basis where 2 n • J is diagonal, 
(IB1I) reduces to a matrix equation, 


(2 jf j+1 


t (2 + 2j, 1 ) 

(2 j - 2f +1 

= — 2 1+2 -' x V \j] 

| t (2 + 2 j, 2 ) 

(~2j + 2) 2j+1 


2 2j— i t (2 + 2j, 2j) 

. (~ij) 2j+1 . 


[ t(2 + 2j,l + 2j) J 


where the Vandermonde matrix for spin j is defined as in (121). So, consider the fcth row on the RHS of (IB2I) : 


2 1+2j f t (2 + 2j, 1) (j + 1 — k) + t (2 + 2j, 2) (j + 1 — k) + t (2 + 2j, 3) (j + 1 — k) 3 + ••• 

(j + l-k) \ +t (2 + 2 j, 2 j) (j + 1 - k ) 2j “ 1 + t (2 + 2 j, 1 + 2 j) (j + 1 - kf j 


If j is semi-integer, say j = n + 1/2 for integer n, then t (2 4 - 2 j, even) = 0, and this fcth row becomes 


2 1+2j 

(j + 1 - k) 


t (2 + 2 j, 1) (j + 1 — k) -f -1 (2 + 2j, 3) (j + 1 — k) 3 + • • • 

+t (2 + 2 j, 2 j - 2) (j + 1 - k) 2 >~ 2 +t( 2 + 2 j, 2 j) (j + 1 - k) 2j 
+t (2 + 2 j, 2j + 2 ) (j + 1 — k) J+ — t (2 + 2j, 2j + 2) (j + 1 — k) J+ 


2 1+2j 

(j + 1 - k) 



t(2 + 2j,2j + 2)(j + l-k) 


2j+2 


(B3) 


(B4) 


where the t (2 + 2 j, 2j + 2) term was added and subtracted to obtain the complete sum on the RHS of (IA3I) . for 
m = j + l/2 = n+l and x = j +1 — k. The sum was then replaced with the product on the LHS of (IA3I) . But 
the product evaluates to zero because one of the terms in the product always vanishes for k > 1 , and therefore 
the &th row on the RHS of (IB2D is 


t (2 + 2 j, 2 + 2j) x (2 j + 2- 2k) 2j+1 = (2 j + 2 - 2 k) 2j+1 , (B5) 

since t (2 + 2j, 2j + 2) = 1. Thus the fcth row on the LHS of (IB2D is obtained, and the identity is established 
for semi-integer j. If j is integer, a similar proof goes through upon using (1A1I) . (See [7], Appendix B.) 
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Appendix C: Biorthogonal matrix examples 

Here are more details about the biorthogonal systems of spin matrices described in the main text, for j = 
1/2, 1, 3/2, and 2. 

Spin j = 1/2 is deceptively simple. The independent powers of the spin matrix are 


S° 


1 0 
0 1 ’ 


s 1 


1 0 
0 -1 ’ 


(Cl) 


and the corresponding trace-ortlronornral dual matrices are the same up to a normalization factor. (NB This is 
not true for any other j.) 


T n = - 


1 

' 1 

0 ' 

1 

' 1 O' 

,, 1 1 

■ 1 1 

2 

0 

1 

’ Tl - 2 

0 -1 

, le - V- =2 

1 -1 


(C2) 


Spin j = 1 is a more interesting example. The independent powers of the spin matrix are given by 



' 1 

0 

0 " 


' 2 

0 

0 


" 4 

0 

0 " 

II 

0 

CO 

0 

1 

0 

, S' = 

0 

0 

0 

, S 2 = 

0 

0 

0 


0 

0 

1 


0 

0 

-2 


0 

0 

4 


(C3) 


and the corresponding trace-orthonormal dual matrices, as well as V 1 , are given by 



' 0 

0 

0 ' 

, T1 = I 

■ 1 

0 

0 

’ — 8 

■ 1 

0 

0 

'N 

II 

00 1 

O 

00 

0 

II 

0 

1 

0 

0 

0 

0 

0 

-2 

0 

2 

0 

-2 


1 

O 

0 

1 

0 

4 

1 

0 

0 

-1 

8 

1 

0 

0 

1 

1 

1 

to 

1 


(C4) 


Spin j = 3/2 is also interesting. The independent powers of the spin matrix are given by 



' 1 

0 

0 

0 " 


' 3 

0 

0 

0 


' 9 

0 

0 

0 ' 


' 27 

0 

0 

0 

S° = 

0 

0 

1 

0 

0 

1 

0 

0 

, s* = 

0 

0 

1 

0 

0 

-1 

0 

0 

, s 2 = 

0 

0 

1 

0 

0 

1 

0 

0 

, s 3 = 

0 

0 

1 

0 

0 

-1 

0 

0 


0 

0 

0 

1 


0 

0 

0 

-3 


0 

0 

0 

9 


0 

0 

0 

-27 


and the corresponding trace-orthonormal dual matrices, as well as V , are given by 


(C5) 


Tn = 



' -1 

0 

0 

0 


■ -1 

0 

0 

0 " 


' 1 

0 

0 

0 " 

1 

0 

9 

0 

0 

1 

0 

27 

0 

0 

m 1 

0 

-1 

0 

0 

16 

0 

0 

9 

0 

II 

001 

0 

0 

-27 

0 

’ T2 "16 

0 

0 

-1 

0 


0 

0 

0 

-1 


0 

0 

0 

1 


0 

0 

0 

1 


(C 6 ) 


T 3 = 



' 1 

0 

0 

0 


' -3 

27 

27 

-3 " 

1 

0 

-3 

0 

0 

, F -1 = — 

-1 

27 

-27 

1 

48 

0 

0 

3 

0 

’ 48 

3 

-3 

-3 

3 


0 

0 

0 

-1 


1 

-3 

3 

-1 


(C7) 


Spin j = 2 helps to establish the general pattern. 


S m = 


4 m 

0 

0 

0 

0 


0 

0 

384 

0 

0 

0 

2 m 

0 

0 

0 

, V- 1 = — 

-16 

128 

0 

-128 

16 

0 

0 

0 

0 

0 

-4 

64 

-120 

64 

-4 

0 

0 

0 

(— 2 ) m 

0 

384 

4 

-8 

0 

8 

-4 

0 

0 

0 

0 

(—4) m . 


1 

-4 

6 

-4 

1 


(C 8 ) 


From the 1st through 5tli rows of V 1 one extracts, respectively, the diagonal dual matrices To through T 4 . 
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Appendix D: A derivation using differential equations 

In this Appendix I derive the coefficients in the Cayley transform spin matrix polynomial by solving differential 
equations. This parallels the first derivation of the CFZ results, as given in [7]. Although the resolvent analysis 
for general matrices, as given in the main text, is much simpler, perhaps the discussion here has some pedagogical 
value in that it emphasizes how much more effort is needed to pursue the same differential equation approach as 
was used for the exponential case. 

For a given spin j let M = 2in • J and write 


1 


2 j 


1 — aM 


= J2 b ™ («) Mm ■ 


Differentiate with respect to a to obtain 


(1 - aM) — 


m —0 


l 


da V 1 ~ ctM 


M 


1 — aM 


That is to say, 


or upon rearranging the terms, 

2 j 


2 j 


db„ 


2 j 


(1 - aM) V ^ M 


m+li 


m =0 


m =0 


2j-l 


E - E Mm+1 ^ (“V) = (“M 


da 


-x± 

da 


(Dl) 


(D2) 


(D3) 


(D4) 


m =0 m =0 

Now apply the fundamental identity dlj to reexpress the RHS of (ID4I) to obtain 

dbo + J2 M m ^ - J2 Mm 4- ( ab m- 1 ) = -J2 Mm (2 i) 2j+1 ~ m t (2 + 2 j, 1+m) A (ab 2j ) . (D5) 


da 


m —1 m—1 m—0 

Trace projections using the dual matrices (i.e. the linear independence of the M m for 0 < m < 2 j) gives the set 
of 1 st order equations 


^ <(2 + 2 i,l )4- ( « M 


db m d 2j+l — m 

*r “ a = - (2!) 


t(2 + 2j, 1 +m) — ( ab 2 j ) for m > 1 . 


(D 6 ) 

(D7) 


These are immediately integrated, with the constants of integration fixed by the requirement that all the b n 
behave properly as a —> 0, namely, bo (a = 0) = 1 and b m (a = 0) = 0 for m > 1. The only additional 
ingredients needed to establish 


bn — 


Trunc [det (1 — 2 ia n • J )] 


(D 8 ) 


det (1 — 2 ia n ■ J) 2 j-r, 

are t (even, odd) = 0 = t (odd, even) [12j 13) and the relations (l39l) and (l40l) . The details are spelled out more 
fully in the following few paragraphs. 

For integer spin (ID61) and (lD7l) become 

dbp 
da 

db 2k 


= 0 


= — (ab 2 k-i) for even m = 2 k , 
da da 

c ^ b ' 2k+1 = (ab 2 k) - 4 J_fc 1 1 (2 + 2 j, 2 + 2k)\ (ab 2 j) for odd m = 2k + 1 , 
da da da 

while for semi-integer spin the equations become 

= —4 j+ 5 \t (2 + 2j, 1)| (ab 2 j) , 
aa aa 


(D9) 

(DIO) 

(Dll) 


db 2k 

da 

db 2 k+i 


= — (ab 2 k~i) - 4 i+ 2 k \t(2 + 2j, 1 + 2k)\ — (ab 2 j) for even m = 2k , 


da 

= (ab 2k ) for odd m = 2k + 1 . 
da 


(D12) 

(D13) 

(D14) 


da 
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Here I have incorporated the phases of the central factorial numbers (see Appendix A) to write 


(D15) 

(D16) 


(— iy k t (2 + 2j, 2 + 2k) = \t (2 + 2 j, 2 + 2k)\ for integer j , 

(i) 2j+1 ~ 2k t(2 + 2j, 1 + 2k) = \t (2 + 2 j, 1 + 2k)\ for semi-integer j . 

Upon integration of (ID 101) and (ID 141) . with the a = 0 initial conditions, it follows that 


b 2 k (a) = ab 2 k-i (a) for 1 < k < j and integer j , 


b 2 k+i (a) = ab 2k (a) for 0 < k < j 


— and semi-integer j . 


Moreover, for integer j 


b 0 (a) = 1 


solves m with the correct value at a = 0. 

Equations (ID11I) and (ID13D also integrate easily to give two difference equations. 

b 2 k+i ~ ab- 2k = -4 J ~ k 1 1 (2 + 2 j, 2 + 2k)\ ab 2j for integer j , 

b 2 k — ab 2 k-i = — U + 2 ~ k 1 1 (2 + 2 j, 1 + 2k)\ ab 2 j for semi-integer j . 


(D17) 

(D18) 

(D19) 


(D20) 

(D21) 


Note the second of these is obtained from the first by k —> k — i. The first of these, along with (ID17I) . may be 
solved by upward recursion starting from k = 0, using (ID19I) . with consistency fixing b 2 j = b 2 j-i. The result 
for integer j is 


bm. — 


a 


ELo 4m « 2m I* ( 2 1 + 2, 2j + 2 - 2m)| 


E 

k =0 


4 k a 2k 


\t(2j + 2,2j + 2-2k)\ 


(D22) 


Note that bo = 1 when j is an integer, and that d m b m /da m | a=0 = ml for all m, a property that is clear from 
the original rational form m- 

Similarly, (ID21I) may be solved by recursion to obtain for semi-integer j, 


bm = —r-r, - 5Z 4 fc a 2fc \t {2j + 2,2j + 2 — 2k)\ . (D23) 

Emir 4 m a 2m 1 1 (2 j + 2 , 2j + 2 - 2 m)| fe=o 

Again note that d m b m /da m | a=0 = ml for all m. But in this case bo ^ 1 since the O ( a 2 -? +1 ) term in the 
denominator is never present in the numerator. Nevertheless, as previously announced [5], in view of (1391) and 
(1401) it follows that the same final determinantal form holds for both integer and semi-integer j, namely (ID 8 I) . 
Technically [33], how sweet it is! 

Alternatively, one may simply substitute (ID22I) into (ID20I) . and (ID23I) into (ID21I) . to verify that those differ¬ 
ence equations are both satisfied, along with (ID 1 71) and (ID18I) . with the proper behavior at a = 0. 

Since for a given j the b n are just renamed ®n J , the result (1341) is proven. 
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Appendix E: Examples of Cayley transforms for j = 1, 2, §, and 3 


det (1 — 2ia n • J )| J=1 / 2 = 1 + < 


1 + 2 ia fi • J 
1 — 2za n • J 


1 + a 2 


((1 - a 2 ) 1 2 x 2 + 4?’a (n • J) 2x2 ) 


det (1 — 2za h • J) | . =J = 1 + 4a 2 


1 + 2za n • J 4za 

1 — 2ia n • J ■ 1 3x3 1 + 4a 2 


((n • J) 

3x3 + 2za (n • J) 3 X 3 ^ 


det (1 — 2za ?7. • «/)|^_g/ 2 — (l “t - a ) (l -7 9a ) — 1 -t- 10a -t- 9a 


1 + 2ia n • J 1 

1 — 2za h • J j _ 3 , 2 (1 + a 2 ) (1 + 9a 2 ) 

—8a 2 

+ (1 + a 2 ) (1 + 9a 2 ) 


((l + 10a 2 — 9a 4 ) I±xA + 4 ia (l + 10a 2 ) (n • J) 4x4 ) 


((n-J) 4 x4 + 2za (n-J) 3 x4 ) 


det (1 — 2za n • J)\^_ 2 = (l + 4a 2 ) (l + 16a 2 ) = 1 + 20a 2 + 64a 4 


1 + 2za fi • J 
1 — 2 ia fi • J J=2 


- *5X5 + (1 + J 2) (1 + ^2) ((" * J )5X5 + (« • J )5X 5 ) 

— 16 za 3 3 4 \ 

a 2 m + 16 a 2 l V (n ‘ J)5x5 + 2la (n ‘ J)sx5 J 


(1 + 4a 2 ) (1 + 16a 2 ) 


det (1 -2ian- J)j J . =5/2 = (l + a 2 ) (l + 9a 2 ) (l + 25a 2 ) = 1 + 35a 2 + 259a 4 + 225a 6 


1 + 2,Oft.J = _1_ M 1 + 35 “ 2 + 259 “ 4 - 225 « 6 > 2 «x« 

1 - 2 ,a ft • J J=5/2 (1 + a 2 ) (1 + 9 a 2 ) (1 + 25 a 2 ) ^ +42q ^ + 35q2 + 259 q 4 } (ft . J) 0x( 

—8a 2 (l + 35a 2 ) / 2 3 \ 

+ fl + a 2 ) (1 + 9a 2 ) FT + 25a 2 ) ( ( ™ ' J)ex6 + (n ' J)ex6 ) 


(1 + a 2 ) (1 + 9a 2 ) (1 + 25a 2 ) 
32a 4 

(1 + a 2 ) (1 + 9a 2 ) (1 + 25a 2 ) 


(i™ ' ^6x6 4" 2za (n • J ) 6x6 ^ 


det (1 - 2za ft • J )| , =3 = (l + 4a 2 ) (l + 16a 2 ) (l + 36a 2 ) = 1 + 56a 2 + 784a 4 + 2304a 6 


1 + 2za- h ■ J I 


= *7x7+ 


4za (l + 56a 2 + 784a 4 ) 


1 — 2za n • J . =3 7X7 (l+4a 2 )(l + 16a 2 )(l + 36a 2 )V v *^7x7^-M 77x7 


n • J) 7x7 + 2za (n • J) 2 


— 16za 3 (l + 56a 2 ) 

(1 + 4a 2 ) (1 + 16a 2 ) (1 + 36a 2 ) 
64za 5 

(1 + 4a 2 ) (1 + 16a 2 ) (1 + 36a 2 ) 


((V- ' **) 7 x 7 4" 2*a (n • J) 7x7 'j 
[(n- J) 7x7 + 2ia (n-J) 6 x7 ) 
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